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Abstract 

We prove that natural generating functions for enumeration of 
branched coverings of the pillowcase orbifold are level 2 quasimod- 
ular forms. This gives a way to compute the volumes of the strata of 
the moduli space of quadratic differentials. 
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^ • 1 Introduction 

■ 1.1 Pillowcase covers and quadratic differentials 
O ! 

lO ■ Consider a complex torus = C/L, where L C C is a lattice. Its quotient 

O 

> ■ 

by the automorphism z —z is a sphere with four (Z/2)-orbifold points 
which is sometimes called the pillowcase orbifold. The map ^ *P is 
essentially the WeierstraB p-function. The quadratic differential (dz)'^ on 
^ ' descends to a quadratic differential on Viewed as a quadratic differential 

■ on the Riemann sphere, {dzY has simple poles at corner points. 
Let be a partition and u a partition of an even number into odd parts. 

We are interested in enumeration of degree 2d maps 

vr : e ^ <P (1) 
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with the following ramification data. Viewed as a map to the sphere, tt has 
profile (z/, 2'^-l'^l/2) over G ^ and profile (2'^) over the other three corners 
of Additionally, tt has profile (/Xj, l^'^^t^i 'j over some £{fi) given points of 
^ and unramified elsewhere. Here is the number of parts in /x. This 
ramification data determines the genus of C by 



In principle, one could allow more general ramifications over and the 
nonorbifold points, but this more general problem is readily reduced to the 

one above^. 

Pulling back (dz)'^ via vr gives a quadratic differential on C with zeros of 
multiplicities {ui — 2} and {2/Xj — 2}. The periods of this differential, by con- 
struction, lie in a translate of a certain lattice. The enumeration of covers tt 
is thus related to lattice point enumeration in the natural strata of the moduli 
space of quadratic differentials. In particular, the d ^ oo asymptotics gives 
the volumes of these strata. These volumes are of considerable interest in 
ergodic theory, in particular in connection with billiards in rational polygons, 
see [6, 18]. Their computation was the main motivation for the present work. 

A different way to compute the volume of the principal stratum was found 
by M. Mirzakhani [19]. 

We thank H. Masur, M. Mirzakhani, N. Nekrasov, R. Pandharipande, and 
A. Zorich for discussions and the anonymous referee for helpful suggestions. 

1.2 Generating functions 
1.2.1 

Two covers tTj : Cj ^ i = 1, 2, are identified if there is an isomorphism 
/ : Ci — > 6*2 such that tti = f 01x2. In particular, associated to every cover vr 
is a finite group Aut(7r). This group is trivial for most connected covers, see 
e.g. Section 3.1 in [7]. We form the following generating function 



^From first principles, the count of the branched coverings does not change if one 
replaces two ramification conditions by the product of the corresponding conjugacy classes 
in the class algebra of the symmetric group. In this way, one can generate complicated 
ramifications from simpler ones. 



x(e) = £(/.)+£(z.)-|/.|-|H/2. 
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where tt ranges over all inequivalent covers (1) with ramification data ji and 
V as above. Note that the degree of any such vr is even. 

In particular, for /j, = v = any connected cover has the form 

: ^ ^ TV± 

with tt' unramified. We have | Aut(7r)| = 2| Aut(7r')| corresponding to the hft 
of ±. Prom the enumeration of possible tt' we obtain, 

n 

By definition, we set 

This enumerates covers without unramified connected components. By the 
usual inclusion-exclusion, one can extract from (3) a generating function 
for connected covers. This generating function for connected covers will be 
denoted by Z°{ii, u; q). 



1.2.2 

Recall the classical level 1 Eisenstein series 

i^2.(g)=^^^^+f:(E^''"'U^ ^=1,2,.... 

n=l y d\n J 

The algebra they generate is called the algebra QM(r(l)) of quasimodular 
forms for r(l) = 5'L2(Z), see [16] and also below in Section 3.3.7. It is 

known that E2, E4, and Eq are free commutative generators of QM(r(l)). 
The algebra QM(r(l)) is naturally graded by weight, where wtE2k = 2k. 
Clearly, for any integer A^, E2k{q^) is a quasimodular form of weight 2k for 
the group 

ro(Ar) = I c = mod ivj C SL2{Z) . 

The quasimodular forms that will appear in this paper will typically be 
inhomogeneous, so instead of weight grading we will only keep track of the 
corresponding filtration. We define the weight of a partition by 

wt /X = + 

The main result of this paper is the following 
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Theorem 1 The series Z'{ji,i>]q) is a polynomial in E'2.{(f), E2{q^), and 
of weight wt II + 11/1/2. 

Several explicit examples of the forms Z'{iJ,, v\ q) are given in the Ap- 
pendix. 



1.2.3 

Quasimodular forms occur in nature, for example, as coefficients of the ex- 
pansion of the odd genus 1 theta-function 

i=l ^ ^ ^ 

at the origin x = 1. The techniques developed below give a certain formula for 
(3) in terms of derivatives of at x = ±1, from which the quasimodularity 
follows. 



1.2.4 

The following discussion closely parallels the corresponding discussion for the 
case of holomorphic differentials in [7, §1.2]. 

Let Q(/x, u) denote the moduli space of pairs (E, 0), where is a quadratic 
differential on a curve E with zeroes of multiplicities {ui — 2, 2fii — 2}. Note 
that we allow i^i = I hence our quadratic differentials can have simple poles. 
For (E, 0) e Q(/x, I/), let E denote the double cover of E on which the differ- 
ential 

cu = 

is well defined. The pair (E,a;) belongs to the corresponding space of holo- 
morphic differentials with zeroes of multiplicity 

{Ui - - - 1} . 

By construction, E is the quotient of E by an involution a. Let P denote the 
set of zeroes of a;; it is clearly stable under a. Then a acts as an involution 
on the relative homology group i/i(E,P, Z). Let H denote the subspace 
of i7i(E,P, Z) on which a acts as multiplication by —1. Choose a basis 
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{7i, ■ ■ ■ , 7n} for H , and consider the period map $ : Q(//, u) — > C" defined 



It is known [18] that 0) is a local coordinate system on Q(/x, i/) and, in 
particular, n = dime = dime Q(Ai, i^). 

PuUing back the Lebesgue measure from C" yields a well-defined measure 
on Q(//, i/). However, this measure is infinite since (j) can be multiplied by 
any complex number. Thus we define Qi(//, v) to be the subset satisfying 



As in the case of holomorphic differentials, the area function is a quadratic 
form in the local coordinates on Q(/i, i/), and thus the image under $ of 
Qi(//, u) can be identified with an open subset of a hyperboloid in C". 

Now let E C Qi(/i, jy) be a set lying in the domain of a coordinate chart, 
and let C^{E) C C" denote the cone over ^(E) with vertex 0. Then, we can 
define a measure p on Qi(/x, i^) via 



where vol is the Lebesgue measure. The proof of [7, Proposition 1.6] shows 
the following analogue: 



where Cov^(/i, u) is the number of inequivalent degree d connected covers 
C — > ^. Thus, the volume p(Qi(/x, i/)) can be read off from the q ^ 1 
asymptotics of the connected generating function Z°{ii, u; q). 

Note that the moduli spaces Q{p, v) may be disconnected. Ergodic theory 
applications require the knowledge of volumes of each connected component. 
Fortunately, connected components of Q(/i, u) have been classified by E. Lan- 
neau [17] and these spaces turn out to be connected except for hyperelliptic 
components (whose volume can be computed separately) and finitely many 
sporadic cases. 



by 





p{E)^Yo\{C<^{E)), 



2D 



p(Qi(/x,z/)) = hm L>-d--°('^''^) J]CovO(/x,i.), 



d=l 
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1.2.5 



The modular transformation 

relates q = and q = 1 and thus gives an easy handle on the g — 1 asymp- 
totics of (3) . This gives an asymptotic enumeration of pillowcase covers and, 
hence, computes the volume of the moduli spaces of quadratic differentials. 

1.2.6 

In spirit. Theorem 1 is parallel to the results of [1, 8, 13], see also [2, 3, 5] for 
earlier results in the physics literature. The main novelty is the occurrence 
of quasimodular forms of higher level. One might speculate whether similar 
lattice point enumeration in the space of A^th order differentials leads to level 
quasimodular forms. Those spaces, however, do not admit an SL2{M.)- 
action and a natural interpretation of their volumes is not known. 

1.2.7 

The following enumerative problem is naturally a building block of the enu- 
merative problem that we consider. Consider branched covers of the sphere 
ramified over 3 points 0,1, oo with profile (z/, 2'^~l'^l/^), (2*^), and respec- 
tively, where /i is an arbitrary partition of 2d. 

The preimage of the segment [0, 1] on the sphere is a graph S on a Rie- 
mann surface (also known as a ribbon graph) with a plenty of 2-valent vertices 
(that can be ignored) and a few odd valent vertices (namely, with valencies 
Ui). The complement of S is a union of £(//) disks (known as cells) with 
perimeters 2^j in the natural metric on S- The asymptotic enumeration of 
such combinatorial objects is, almost by definition, given by integrals of ip- 
classes against Kontsevitch's combinatorial cycles in Mg^£(^), see [15]. There 
is a useful expression for these integrals in terms of Schur Q-functions ob- 
tained in [4, 11]. In fact, our original approach to the results presented in 
this paper was based on these ideas. 

While the proof that we give here is more direct, it is still interesting to 
investigate the connection with combinatorial classes further, especially since 
a natural geometric interpretation of combinatorial classes is still missing. 
Perhaps the Gromov-Witten theory of the orbifold ^ is the natural place to 
look for it. This will be further discussed in [22]. 
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2 Character sums 



2.1 Characters of near- involutions 



2.1.1 



There is a classical way to enumerate branched coverings in terms of irre- 
ducible characters which is reviewed, for example, in [10] or in [21]. Special- 
ized to our case, it gives 



Z{^^,.■q) = Y,q\^\'- ^ f.2,2,...(A) f,,,,...(Af J] ^ (4) 



where summation is over all partitions, dim A is the dimension of the cor- 
responding representation of the symmetric group, and fr;(A) is the central 
character of an element with cycle type rj in representation A. Recall that 
the sum of all permutations with cycle type rj acts as a scalar operator in any 
representation A and, by definition, this number is f^(A). In (4), as usual, we 
abbreviate fk,i,i,... to f^. 



A lot is known about the characters of the symmetric group S{2d) in the 
sitTiation when the representation is arbitrary but the support of the permu- 
tation is bounded by some number independent of d. In particular, explicit 
formulas exist for the functions f^. 

Understanding the function iy^2,2,... is the key to evaluation of (4). That 
is, we must study characters of permutations that are a product of a per- 
mutation with finite support and a fixed-point-free involution. We call such 
permutations near-involutions. 



By a result of Kerov and Olshanski [14] , the functions fjt belong to the algebra 
A* generated by 



p,(A) = (1 - 2-^)a-k) + ^ [(A, -i + \f- i-i + m , (5) 



moreover has weight k + 1 in the weight filtration on A* defined by setting 

Wt Pfc = /c -|- 1 . 



A 




2.1.2 



2.1.3 
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The functions are central characters of certain distinguished elements in 
the group algebra of symmetric group known as completed cycles. See [21] 
for the discussion of the relation between and from the viewpoint of 
Gromov-Witten theory. 



Our next goal is to generalize the results of [14] to characters of near- 
involutions. This will require enlarging the algebra of functions. In addition 
to the polynomials p^, we will need quasi-polynomial functions p^ defined in 
(6) below. 

It is convenient to work with the following generating function 



2.1.4 



e(A, z) J2 e^('^-^+^) 




k 



By definition, set 



Pfe(A) = i k\ [z^] e(A, z + m) 

= E - ^ + i)' - + \f] + const , 



(6) 



where the constant terms are determined by the expansion 




k 



Up to powers of 2, they are Euler numbers. 



2.1.5 



Define 



A = Q [Pik, Pfe] 




where A' denotes the conjugate partition. 
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2.1.6 



In the definition of A, we excluded the function 

i 

which measures the difference between number of even and odd numbers 
among {Xi — i + 1} , also known as the 2- charge of a partition A. 

Every partition A uniquely defines two partitions a and /3, known as its 
2-quotients, such that 

{A,-i + i} = {2(Q;,-i + i)+po(A)} U {2(A-i + i)-po(A)} . 

A partition A will be called balanced if po(A) = |. 

Several constructions related to 2-quotients will play an important role 

in this paper. A modern review of these ideas can be found, for example, in 
[9]. In particular, it is known that the character X22 ^ fixed-point free 
involution in the representation A vanishes unless A is balanced, in which 
case 

) dima dim/3 . (8) 
\P\J 

It follows that only balanced partitions contribute to the sum (4) . 




2.1.7 

For a balanced partition A, define 

g„(A) = 5?«4^. (9) 

H2,2,...)[X) 

We will prove that this function lies in A in the following sense: 

Theorem 2 The ratio (9) is the restriction of a unique function g,^ e A 0/ 
weight \v\l1 to the set of balanced partitions. 

Several examples of the polynomials gu can be found in the Appendix. 
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2.1.8 



In view of Theorem 2, it is natural to introduce the following pillowcase 
weight 

w(A)=(^|^)%«JAr. 
Theorem 1 follows from (4), Theorem 2, and the following result: 

Theorem 3 For any F E A, the following average 

A 

is a polynomial in E2{q^), £'2(5^), and Ei{q^) of weight wt F . 

Note that if F is homogeneous of odd weight then {F)^ = 0. This can 
be seen directly from (7). Also note that (10) will not in general be of pure 
weight even if F is a monomial in the generators and pfe. This contrast 
with [1, 8] hints to the existence of a better set of generators of the algebra 
A. Probably such generators are related to descendents of orbifold points in 
the Gromov-Witten theory of 

2.1.9 

It will be convenient to work with the following generating functions for the 
sums (10) 

F{x^,...,Xn) = (j{e{\Mxi))^. (11) 
The function (11) will be called the n-point function. 

2.2 Proof of Theorem 2 
2.2.1 

In the proof of theorems 2 and 3 it will be very convenient to use fermionic 
Fock space formahsm. This formahsm is standard and [12, 20] can be rec- 
ommended as a reference. A quick review of these techniques can be found, 
for example, in Section 2 of [21]. We follow the notation of [21]. 
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2.2.2 



By definition, the space qV is spanned by the following infinite wedge 
products 

- Ai - I A A2 - I A A3 - I A ■ ■ ■ , (12) 

where ^, /c e Z + | is a basis of the underlying space V and A is a partition. 
The subscript in A^q^ refers to the charge zero condition: the ith factor 
in (12) is — i + I for all sufficiently large i. 

There is a natural projective representation of the Lie algebra gl(V^) on 
A^o^- For us, the following elements of qI{V) will be especially important: 

ek[f(x)]i^f(i-l)i^, (13) 

where / is a function on the real line. To define the action of £0 [/(a^)] on 
A^oV one needs to regularize the infinite sum J2i<ofil ~ 0- This regular- 
ization is the source of the central extension in the qI{V) action. When / is 
an exponential as in 

£,k{z) = Sfc [e^l , 

this infinite sum is a geometric series and thus has a natural regularization. 
By differentiation, this leads to the ^-regularization for operators £.k[f] with 
a polynomial function /. 

2.2.3 

Other very useful operators are 

a;ifc = £fe[l], ky^O. 

The operator H defined by 

Hvx^ \X\vx 

is known as the energy operator. It differs only by a constant from the 
operator S,o[x]. The operator H defines a natural grading on A'^'qV and 
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2.2.4 

A function F{X) on partitions of n can be viewed as a vector 

\X\=n 

of energy n. For example, the vectors 

i/^) = ^n«---0=^E^>A, (14) 

where 

M = I Aut/x| JJ/Xj , 

correspond to irreducible characters normalized by the order of the central- 
izer. 

2.2.5 

The operator E.o{z) is the generating function 

z '-^ k\ 

k 

for the operators CPfc acting by 

In parallel to (6), we define operators 5*^ by 



:£o(^ + 7ri) = ^y'yk- 



k 



Translated into the operator language, the the statement of Theorem 2 is 
the following: the orthogonal projection of |z/, 2'^~l^l/^) onto the subspace 
spanned by the vx with A balanced is a linear combination of vectors 

W^..]!^-..]^') (15) 

with 

wt ^ + \j2\ < \v\l2 
and coefficients independent of d. 
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2.2.6 

Let us call the span of vx with A balanced the balanced subspace of A^o^- 
A convenient orthogonal basis of it is provided by the following vectors 

|p; p) = n ^-pi n °^-pi ' Pi,Pie2Z, (16) 

where the operators ak are defined by 

a, = i'+'Ek{7Tt) = J](-l)"+^i^„_,,„ + I , (17) 

n 

the operators E^j being the matrix units of qI{V). From the commutation 
relations for the operators Ek{z) we compute 

[ak, arr^] = [(-1)*= - (-I)'"] + k{-l)%+m , (18) 

K, a^] = [1 - {-!)'] (^ak+m + . (19) 

In particular, when both k and m are even, all these operators commute 
apart from the central term in [oik, ci;_fe]. 
The adjoint of cifc is 

al = (-l)'=a-fc, 
which gives the following inner products 

provided all parts of all partitions in (20) are even. In particular, the vectors 
(16) are orthogonal. It is clear that they lie in the balanced subspace and 
their number equals the dimension of the space. Therefore, they form a basis. 

2.2.7 

The projection of |z/, 2'^^l'^l/^) onto the balanced subspace is given in term of 
inner products of the form 

(l/,2'^-|''l/2|(p,2'^-|''l/2-|^l/2);p> 
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where all parts of u are odd, all parts of p and p arc even, and p has no parts 
equal to 2. From the commutation relations (18) and (19) we conclude that 
this inner product vanishes unless 

The nonvanishing inner products are the following 

2t(y)-m 

where the combinatorial coefficient C(i/, p) equals the number of ways to 
represent the parts of p as sums of parts of v. For example 

C((3,l,l,l),(4,2)) = 3, C((3,l,l,l),(6)) = l. 

2.2.8 

The matrix elements 

i2'|n^M^n^/^.|('^'2'"''''^'"'''^^ -^^^2, (22) 

describe the decomposition of the vectors (15) in the basis (16). Since 

3^1 |2'^> = (d-^) |2<^> , (23) 

we can also assume that pi ^ 1. 

We claim that (22) vanishes unless 

wt/x+ > wtp/2+ |p|/2, (24) 

where p/2 is the partition with parts pi/2 (recall that all parts of p are even). 

2.2.9 

The usual way to evaluate a matrix element like (22) is to use commutation 
relations to commute all lowering operators to the right until they reach the 
vacuum (which they annihilate) and, similarly, commute the raising operators 
to the left. 

We will exploit the following property of the operators (Pfe and 7k'- their 
commutator with enough operators of the form a^2pi and Q;_2pi vanishes. All 
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such commutators have the form £fe[/] with f{x) = (±l)^p(x), where p{x) 
is a polynomial. Commutation with a-2p, takes a finite difference of p{x); 
commutation with a-2pj additionally fiips the sign of ±1. 

Since a (/c + l)-fold finite difference of a degree k polynomial vanishes, 
the commutator of of Tfe with more than k + 1 operators of the form a-2pi 
or Q;_2pi vanishes. In fact, a (A; + l)-fold commutator may be nonvanishing 
only because of the central extension term. To pick up this central term, the 
total energy of all operators involved should be zero and the number of a's 
should be even. Same reasoning applies to T^, but now the the number of 
a's should be odd to produce a nontrivial {k + l)-fold commutator. 

2.2.10 

Now look at one of the raising operators involved in (22), say a_p.. This 
operator commutes with eta and its adjoint annihilates the vacuum, so only 
the terms involving the commutator of a-p^ with one of the or give a 
nonzero contribution to (22). The commutator [J'^.,a_p.] has energy {—pi) 
and so its adjoint again annihilates the vacuum. The same is true for the 
commutation with 7p^ . To bring these commutators back to zero energy, one 
needs to commute it pi/2 times with 02 • Given the above bounds on how 
many commutators we can afford, this implies (24). 

2.2.11 

When the bound (24) is saturated, then a further condition 

£{p)+£{p)>£{l,)+£{fl) 

is clearly necessary for nonvanishing of (22). The unique nonzero coefficient 
saturating both bounds corresponds to 

P = 2// , p^2j2. 

Moreover, when divided by the norm squared of the vector | (p, 2<^-I''I/2-IpI/2) ; , 
this coefficient is independent of d. 

2.2.12 

For general p and p, the similarly normalized coefficient will be a polynomial 
in d of degree 

^(wt/i+|/x| -wtp/2- |p|/2) (25) 
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because so many operators q;_2 can commute with "P^^s or ^pi^s instead of 
commuting directly with 02 's. 

By induction on weight and length, we can express the basis vectors (16) 
in terms of (15) with 1 and coefficients being polynomial in d of degree 
at most minus the difference (25). By (23), to have d-dependent coefficients 
and 7^ 1 is the same as to allow ^ii — l and make coefficients independent 
of d. The bound of degree in d ensures that this transition preserves weight. 
This concludes the proof of Theorem 2. 

3 Proof of Theorem 3 

3.1 The pillowcase operator 
3.1.1 

Consider the following operator 




Because this operator is normally ordered, its matrix elements (SBv, w) are 
well defined for any vectors v and w of finite energy. The relevance of this 
operator for our purposes lies in the following 

Theorem 4 The diagonal matrix elements of SB are the following: 

f w(A) , A is balanced 

(SB'Ua,va) = { 

I , otherwise . 

The proof of this theorem will occupy the rest of Section 3.1. 
3.1.2 

Let N be chosen so large that A2Ar+i = 0. Then because the operator (26) 
is a product of an upper unitriangular and lower unitriangular operator, the 
vectors Aj — i + | with i > 2N in 

^^A = Ai - I A A2 - I A As - I A . . . 
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are inert bystanders for the evaluation of (SBi'a, vx)- The whole computation 

is therefore a computation of a matrix element of an operator in a finite 
exterior power of a finite dimensional vector space V^^"^ with basis 

ek^-2N + k + l, k^0...,Xi + 2N -1. 

By definition, matrix elements of SB in exterior powers of yl^l are determi- 
nants of the matrix elements of SB acting on the space V^I^l itself. The latter 
matrix elements are determined in the following 



Proposition 1 We have 
h{k)h{l) 

where 



1 , k = l = mod 2 

, k = l = l mod 2 (27) 

2/{k — I) , otherwise , 



b(^) 



2k lk/2\\ 



3.1.3 

For the proof of Proposition 1, form the following generating function 

k,l 

From the equality 
and definitions, we compute 



\ / i + X 



2n + 1 / \ 1-x 



1 — xy V 1 — X V 1 + y 
The following factorization 

d d\ {x + y)il + x)il-y) 



dx "dyj-'^'^' (1 -x2)3/2(l_ ^2)3/2 ' 
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by elementary binomial coefficient manipulations proves (27) for k ^ I. To 
compute the diagonal matrix elements observe that the above differential 
equation uniquely determines f{x,y) from its values on the diagonal x = y. 
On the diagonal, the skew-symmetric terms in (27) cancel out and evaluation 
is immediate. 



We now proceed to the computation of the matrix element {^vx,v\). We 
have the following 

Proposition 2 We have 



provided A is balanced and {^v\, vx) — otherwise. 

The proof of this proposition is the following. Observe that by Proposition 
1 the matrix element (SBt'A, ^a) is a determinant of a 2N x 2A'" block matrix 
in which the odd-odd block is identically zero, the even-even block is a rank 1 
matrix with all elements equal to 1 and the off-diagonal blocks have the form 



y ^^y j , where {xi] and {yi] are the odd and even subsets of {Aj — i + 2N}. 
Since the odd-odd block is identically zero, its size has to be < for the 
determinant to be nonvanishing. Similarly, if the size of the even-even block 
is larger than N then the determinant is easily seen to vanish. It follows 
that both blocks have size N , which precisely means that the partition A is 
balanced. It remains to use Cauchy determinant 



Note that decomposition of {\i — i + 2N} into the even and odd subsets is the 
same as the 2-quotient construction from Section 2.1.6. Theorem 4 follows 



3.1.4 






to finish the proof. 



3.1.5 
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from formula (8) and the following classical formula 

dim A _ Ili<j<Ni^i - + J - '>') 
|A|! " ll[X^ + N-t)\ 

where N is any number such that Ajv+i = 0. 
3.1.6 

It would be interesting to find an interpretation of the operator SB in confor- 
mal field theory. Note that 

\n>0 / \ n>0 / 

is the Wiener-Hopf factorization of the function taking the value on the 
upper /lower half-plane. 

3.2 Formula for the n-point function 
3.2.1 

Theorem 4 yields the following operator formula for the n-point function (11) 

F(xi , . . . , x„) = z{^^-q) n ® ' (28) 

where the trace is taken in the charge zero subspace of the infinite wedge and 
H is the energy operator 

We have the following expression for the operator £o in terms of the fermionic 
currents: 

£o(lna;) = [y'^]il){xy)il)*{y) , 

where [y^] denotes the constant coefficient in the Laurent series expansion in 
the variable y. Therefore, 

F(xi,...,x„) = ^^^^^_^^x 

[yl ... 2/°] tr g^^(a;i2/i) rivi) ■ ■ ■ H^nVn) riVn) 3B . (29) 
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3.2.2 

By the main result of [9] , we have 

w(A) < 1 (30) 

for any partition A. In other words, all diagonal matrix elements of SB are 
bounded by 1. For the off-diagonal elements, we prove the following cruder 
bound 

Proposition 3 Let M — max{|A|, |//|}. Then 

(3B vx, v^) < exp I - J2 I ^ const M^^ . (31) 

To see this note that 

where SB'^^ is the following truncated operator 



exp 



K2n+1<M / \ 2n+l<M 



We claim that the operator SB^^^' is a multiple of a unitary operator. Indeed, 

+ 



\ n=0 / 



whence the result. 

In fact, we will only use that (31) is bounded by a polynomial in the sizes 
of the partitions. 

3.2.3 

By normally ordering all fermionic operators in (29) and using the estimate 
(31) one sees that the trace converges if 

\yn/q\ > \xiyi\ > \yi\ > ■ ■ ■> IXnVnl > IVnl > 1 ■ (32) 
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3.2.4 

The proof of the following identity is given in [12], theorem 14.10: 



/ v-^ (xy)"' — I/" \ / v-^ y — ixy) " \ 
exp I «_„ I exp I 2^ an | . (33) 



It allows to express the operator in (29) in terms of bosonic operators 

With respect to the action of the operators a„, the charge zero subspace of 
the infinite wedge space decomposes as the following infinite tensor product 



oo oo 



n=l fe=0 



the distinguished vector in each factor being v^. This gives a factorization of 
the trace in (29). The trace in each tensor factor is computed as follows: 



^fc 1 — q 



1 fnABq'' 
exp 



For example, this shows that 

trg^38=(g2)-V2 = Z(0,0;g), 

where 



/"nAB 
V 1 - J 



(a)oo = n(l-0' 



n>0 



and so the 0-point function is F{) = 1, as expected. For the n-point function 
this gives the following 



Theorem 5 We have 



F{x,,...,Xn) = l[-^x 



where the series expansion is performed in the domain (32) . 
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3.3 Quasimodular forms 
3.3.1 

In the computation of (34), we can assume that 1 < \xi\ -C \q~^\ for all i and 
hence 

The series expansion in (34) can then be performed using the following ele- 
mentary 

Lemma 4 We have 

provided c> \bi\ > \q\c for i = 1, . . . ,n. 

This is obtained by computing the difference of ^j^i^^, and ^^^^^^^^ as a 
sum of residues using 

^'{1) = 1. 

3.3.2 

There are two obstacles to literary applying this lemma to the evaluation 
of (34). The first are the square roots in (34). However, we are ultimately 
interested in the expansion of (34) about Xi — ±1. The expansion of the 
integrand about Xi — ±1 contains no square roots, only theta function and 
its derivatives. Formulas for integrating derivatives can be obtained from 
(35) by differentiating with respect to parameters. 



3.3.3 

The other issues is that at = 1 the integrand is an elliptic function of the 
corresponding yj, and so the left hand side of (35) gives infinity times zero. 
This can be circumvented, for example, by replacing each factor of Xi in the 
argument of each theta function an independent variable and specializing 
them all back to Xi only after integration. By the L'Hospital rule, this will 
produce an additional differentiation any time we expand around Xi = 1 for 
some i. 
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3.3.4 



In the end, we will get some rather complicated polynomial in theta functions 
and their derivatives evaluated at ±1 divided by a power of !)■ This 
means that we will get a combination of Eisenstein series arising from 

2k 

(36) 

and 

^ = 2 E (J [E'^M - -i'^E.M')] , (37) 
together with the following product 

Note that (38) has weight —1. 
3.3.5 

Without knowing the precise form of the answer, one can still make some 
qualitative observations about it. 

Suppose we are interested in the coefficient of . . . z^" in the expansion 

of 

in powers of Zi. We claim that the weight of this coefficient is at most Yl ki+r. 
Indeed, we from (36) and (37) we have 

wt I X— I 'ii}(x) = k — 1 . 
\ ax J x=±i 

This gives the following count for the weight 

n-n + ^ ki + r , 

where the first n is added because of the prefactor in (34), the second n is 
subtracted due to integration in Ui (which, by Lemma 4 changes the balance 
of ^-factors by 1), is the number of times we need to differentiate the 
integrand, and, finally, r additional differentiations are needed for reasons 
explained in Section 3.3.3. 
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3.3.6 



We further claim that (34) is, in fact, a polynomial in the coefficients of (36), 
(37) and 

^ = - ''EM - + l^EM) . (39) 

First observe only even powers of (38) appear in the answer. This is because 
the formula (34) has a balance of minus signs in the arguments of theta 
functions in the numerator and denominator. Every time we specialize t/j 
to one of the poles in (35), the balance of minus signs changes by an even 
number. 

Inverse powers of (39) cannot appear in the answer because they grow 
exponentially as g ^ 1 and there arc no other exponentially large terms to 
cancel this growth out. The averages (10) may grow only polynomially as 
g — > 1 because of the bound (30). 



3.3.7 

Recall from [16] that a quasimodular form for a congruence subgroup F C 
SL2{'L) is, by definition, the holomorphic part of an almost holomorphic 
modular form for F. A function of |g| < 1 is called almost holomorphic if it 
is a polynomial in (ln|g|)~^ with coefficients in holomorphic functions of q. 
Quasimodular forms for F forms a graded algebra denoted by QM(F). By a 
theorem of Kaneko and Don Zagier [16], 

QM(F) = Q[^2] ® M(F) . 

In particular, 

EM, EM), E2{q^) e QM(Fo(4)) (40) 

where 

ro{A)^^(l c = 0mod4| C ,SL2(Z). 

Hence all averages (10) lie in QM(Fo(4)). 

In fact, the series (40) generate the subalgebra QM2*(Fo(4)) of even weight 
quasimodular forms. This is because M2*(Fo(4)) is freely generated by two 
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generators of weight two, for example, by E2^^{q) and E2^^{q'^), where 
Er{q)=E,{q)-2E,{q') = ^ + f2( ^J^'^- 

n=l \d|n,dodd / 

3.3.8 

Note that because w(A) = for any partition A of odd size, the series (10) 
is in fact a series in q'^. It follows that it is quasimodular with respect to a 
bigger group, namely 

(J °)r„(2) r,4). 

In other words, (10) is, in fact, obtained by substituting q ^ (f into an 
element of QM(ro(2)). We have 

M{T^{2)) = Q[Ef\q),E^{q^)] 

and hence 

QM(ro(2)) = Q[E2{q), E^iq^), E^{q^)] . 
This concludes the proof of Theorem 3. 

A Examples 

In this appendix, we list some simple examples of the quasimodular forms 
Z'(/x, I/; q) appearing in Theorem 1 and polynomials gj, from Theorem 2. 

A.l Quasimodular forms Z'{ii^ v; q) 

Z'((l, 1), (2)) = 20£;2(g^)' - 2QE2{q') E^iq') ^AE^iq^f - 5/3E,{q^) . 



^'((3,1), (3)) 



2112 „ , ^,3 3888 „ , /i.2„ / Ox 2304 „ , a. „ , 0^2 
E2{qy + ^ E2{q^) E^iq^) — E^iq") E2{q^) 



+ 



384 



E2{q^f + ASE^q^) ^2(9^) - 36^4(9^) ^2(9' 
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Z'((3, 3), (2)) = E,{q'f - E,{q'fE,{q^) 

+ ^E2{q')E2{q'f E,{q'f - 2AE,{q') E,{q') + ?> E,{q') E,{q') 
5 5 

+ 1^/2 E^iq'f - 1^/2 E^iq") E^iq^) + 3/2E2{q^f - 5/8E^{q^) 



Z'((5, 1), (2)) = ^ - 1160 E,(q'fE,{q') + 280 £;2(g^)' 

- f E,{q^f - ^ E,{q') E,{q') + f E,{q') E,{q^) + ^-f E,{q'f 

- ^ E,{q') E,{q^) + | E,{q^f - ^ E,{q') . 

Z'((l,l,l,l),0) = l/4E2(g^) + l. 



Z'((3, 3, 1, 1), 0) = ^ - 12 + y i^2(g') i?2(g') 

- 9/4 + 5/4 £;4(g^) + ^ E,{q') + 3/8 E^iq') . 



Z'iiB, 1, 1, 1), 0) = ^ - 10 E,iq'f + 15 

-5/2E,{q'f+^^E,{q') + ^E,{q'). 



Z'((3,3,3,3),0) = 



39 



E,{q^f - ^ E,{q'fE,{q^) + ^ £;2(g^) ^J^l?^)' 
zO 



7 



33 



141 



10 



16 

32 ^ 32 ^ 256 ^ ^ 32 ^ 2048 
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Z'((5, 3, 3, 1), 0) = 132 EM'f - 708 E^'f E^{q^) + 639 E^{q') E,{q'f 

- 114 E2{q^f - 15 E^{q^) E^iq*) + 55 E^{q^) E^iq^) - 310 E2{q''f 
1365 ^ , A. ^ , 285 „ , 9,2 175 „ , 615 „ , 4, 85 „ , 9, 375 

A. 2 Polynomials gjy 

gi,i = l/2pi. 
g3,i = l/6p2 + l/6p2-l/2pi. 

g3,3 = P? + 1/18 P1P2 + ^ P3 - 1/4 P2 + 3/16 pi . 



25 

g5,i = 1/30 p? + 1/10 P1P2 - 1/2 pipi + 1/15 p3 - 1/2 P2 + ^ Pi 



1/1 ^/^^o 2 25 



55,3 



Pi - P?P2 - 1/12 p^pi + — P3P1 + 7;^ Pi + 1/40 pI 



360 "' 60 "^'^ ' 45"^ 36 

9 1 , 25 25 

- 1/12 P2P1 + 5/8 pi^ + _ p4 - 1/2 Pa + — p2 - — pi 



gi,i,i,i = -l/24p? + l/12p2 + ^. 



g3,l,l,l = Pi ~ 1/2^ P1P2 ~ 1/4 PlPl + ^ P3 + 2/^ Pi 



1 1 9 1 

g3,3,l,l = ^ Pi ~ 1/12 PlPl + ^ P3Pl ~ 1/8 P2P1 + ^ Pi ~ 7^ P2 

- 1/12 P2P1 + 1/8 pi' + 1/36 P4 + 4 P2 - 3/4 Pi + 7^ 

16 256 
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g3,3,3,l 
1 



4860 



Pl + 



486 



P?P2 + 



108 



P?Pi 



972 



PaPi 



l/24p2P?--p? 



+ P1P2 - 1/36 P1P2P1 + — P4P1 - 7^ P3P2 - — P3P1 - 1/24 P2P2 



324 



972 



108 



162 

31 _ _ 19 _ _ 2 _ 153 

+ 1/8 P2P1 + — P1P2 - — PlPl + 1/4 P3 - P2 + 77^ P5 + 777^ Pi ■ 



96 



32 



405 



128 



g3,3,3,3 - ^ ^'-2^^'^' + ^^'^'-^^'+!^^'^'^-^ P^Pl 
+ ^ P3P1P2 - ^ P2P1P2 - ^ P?P2 - 1/12 P?P, - ^ p3 - -1- p^ 

1 - 59 _ _ o 3 _ 1 _ _ 231 , _^ 

' 216 P^P^ + 864 P^P^ + P^ - 64 P^P^ + 1215 P^P^ + 512 P^ + P^ 
2 5 _ 5 1 _ 129_ 9 27 



1/12 P2P1 + 3/8 pi + — P4 - 777 P3 + 777TT^ P6 + 77^ P2 - o Pi + 



144 12^"^ 2916 '^"^ 256" 8" 2048 
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